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NOMENCLATURE 


The  definitions  of  the  following  symbol*  are  used  throughout  the 
report  except  where  otherwise  stated.  Some  symbols  defined  in  the 
text  of  the  report  are  not  included. 

* 

a . Velocity  of  wave  propagation  in  a  liquid- filled  pipe.' 


A 

*  *  •  •  •  •  • 

b . 

8  •  •  •  •  • 

B-Cju)  .  .  . 

.  .  . 


Cross  sectional  area  of  pipe. 

Thickness  of  pipe  wall. 

Constant  used  in  determining  the  frictional  resistance 
to  flow. 

System  function  for  the  ith  pipe  at  point  x^  . 

System  function  for  the  ith  pipe  at  Xj  =  i.  . 

Velocity  of  wave  propagation  in  an  infinite  medium. 
Complex  Fourier  coefficient. 


C 

D 


e 

E 


(Xj.  **>)'• 


Eqk(xi.M. 


f 


Coefficient  of  capacitance. 

Inside  diameter  of  pipe. 

Base  of  natural  logarithms  =  2.  718*  •  •  . 

Modulus  of  elasticity  of  pipe  material. 

Transfer  function  relating  pressure  at  point  x.  to 
pressure  at  point  k  . 

Transfer  function  relating  pressure  at  point  xi  to  flow 
rate  at  point  k  . 

Fundamental  frequency. 

Cut-off  frequency. 


JDarcy-Weisbach  friction  factor. 
[Cyclic  frequency. 


iii 


Discrete  frequency. 


F  .  (x.,jw) .  Transfer  function  relating  power  at  point  x.  to  pressure 
pK  1  at  point  k  .  1 


F  ,(x.,jw).  Transfer  function  relating  power  at  point  x. 
qK  1  at  point  k  .  1 


to  flow  rate 


G.(jcu) 


G(fj  . 
h  (t)  . 
H^) 


HJ(jw) 


Jq  and 


qi  * 


K,i  • 


.  Transfer  function  relating  pressiyre  at  point  x.  to  flow 
rate  at  terminus.  1 

.  Transfer  function  relating  pressure  at  point  i.  to  flow 
rate  at  terminus. 

+ 

.  Spectral  density  as  defined  by  American  Standard;-  - 

% 

.  System  weighting  function. 

.  Transfer  function  relating  flow  rate  at  point  x.  to  flow 
rate  at  terminus. 

.  Transfer  function  relating  flow  rate  at  point  i.  to  flow 
rate  at  terminus. 


(Refers  to  ith  segment  of  pipe. 
General  subscript. 


Bessel  functions  of  the  first  kind. 


General  subscript. 

Bulk  modulus  of  a  liquid. 

Gain  function  relating  pressure  at  point  x.  to  flow  rate 
at  point  xA  =  0  .  1  ^ 

Gain  function  relating  flow  rate  at  point  x.  to  flow  rate 
at  point  =  0  . 

Gain  function  relating  flow  rate  at  point  i.  to  flow  rate 
at  point  x  .  «  0  . 

Gain  function  relating  pressure  at  point  x.  to  flow  rate 
at  terminus. 

Gain  function  relating  flow  rate  at  point  x.  to  flow  rate 
at  terminus . 

Length. 

Length  of  ith  pipe  segment. 


iv 


p  ft 


4 


4 


* 


L 


m 


M.(jco)  .  .  . 


Coefficient  of  inertia.. 

One-sided  Laplace  transform. 

Two-sided  Laplace  transform. 

r~  ' 

Number  of  pipes  in  series.  . 

Number  of  pipes  in  parallel. 

General  subscript. 

System  function  for  the  ith  pipe  at  point  . 


Mj(jw)  .  .  . 


n 


P 


System  function  for  the  ith  pipe  at  x^  =  . 

_  .  *  JatV* 

[Exponent  of  q  for  mean  flow. 

<Refers  to  source  en<$. 

[General  subscript. 

Pressure  variation  as  a  function  of  distance  and  time. 


P 


Mean  pressure  as  a  function  of  distance. 


Pt . 

P(X.,j(0)  .  . 

•  •  • 
fc(X,  s) .  .  . 


Instantaneous  pressure  as  a  function  of  distance  and 
time  =  p  +  p  . 

Pressure  variation  at  point  x.  transformed  with  respect 
to  time  with  s  s  jw  . 

Transformed  pressure  variation  at  x^  =  0  with  s  =  ju>  . 

Pressure  variation  transformed  with  respect  to  time  and 

di  stance . 


q 


Volume  flow  rate  variation  as  a  function  of  distance  and 
time. 


q .  Constant  mean  volume  flow  rate. 


qt . Instantaneous  volume  flow  rate  as  a  function  of  distance 

and  time  ”  q  +  q  . 

Q(x.,jo)  .  .  Volume  flow  rate  variation  at  point  x^  transformed  with 
1  respect  to  time  with  s  =  jw  . 


Qj(jw)  .  .  .  Transformed  flow  rate  variation  at  x^  =  0  with  s  =  ju>  . 


B)  • 


Flow  rate  variation  transformed  with  respect  to  time 
and  distance. 


r 


Refers  to  terminus. 


R 


Coefficient  of  resistance. 


v 


s . Laplace  transform  with  respect  to  time  operator. 

t . Time. 


T  jLimit  of  a  period. 

. ^Denotes  truncation. 

Tjfju)  .  .  .  System  function  for  the  ith  pipe  at  point  xi  . 

Tj'(j<*>)  •  •  •  System  function  for  the  ith  pipe  "it  x^  -  i.  . 

u. . Instantaneous  fluid  velocity  as  a  function  of  distance  and 

time. 

m 

w  .....  .  Power.  ..  *..**•'■•* 

*  ; 

W(Xj,jw).  .  Transformed  power  *at  point  x^  with  s  -  jw  . 

Wave*  •  •  •  Average  power, 
x  .....  Axial  distance. 

Xj . Axial  distance  measured  from  the  terminating  end  of  the 

ith  segment. 

x(t)  ....  Function  of  time. 

X(s)  ....  Laplace  transform  of  x(t)  . 

Y  .(xi»  jw)*  Transfer  function  relating  flow  rate  at  point  x.  to 
P  pressure  at  point  k  . 

z . Characteristic  value. 

Z(x.,  j<*>)  .  .  Ratio  of  transformed  pressure  to  flow  rate  at  point  x^  . 

Zj(jw)  .  .  .  Ratio  of  transformed  pressure  to  flow  rate  at  x^  =  0  . 

ZQ . Characteristic  impedance. 

zci  •  •  •  •  Characteristic  impedance  of  the  ith  segment. 

Zr . Ratio  of  transformed  pressure  to  flow  rate  at  terminus. 

Z  .  (x.,  jw) .  Transfer  function  relating  pressure  at  point  x.  to  flow 
q  rate  at  point  k  . 

or 
*i 


Attenuation  constant. 

Attenuation  constant  of  the  ith  segment. 


0 . Phase  constant. 

0j . Phase  constant  of  the  ith  segment. 


Propagation  coefficient. 

/ 

Propagation  coefficient  of  the  ith  segment. 

Reflection  coefficient. 

Reflection  coefficient  of  the  ith  segment. 


0 . Phase  angle  of  reflection  coefficient. 

X . Laplace  transform  with  respect  to  distanc'6  operator! 

* 

p . Mass  density. 

(Dummy  time  variable. 

. ^Displacement  in  tijpe. 


Phase  angle  of  the  ratio  of  pressure  to  flow  at  a  point  k 


0  (t)  .  .  .  Autocorrelation  function  of  x(t)  . 

AA 


0  (t)  .  .  .  Cross-correlation  function  of  x(t)  and  y(t)  . 

xy 

0pq(T)  *  •  *  8011,111  power. 

^xx^01*)  •  Spectral  density  of  x(t)  . 

$  (ju)  .  .  Cross-spectral  density  of  x(t)  and  y(t)  . 
xy 

$  (jw)  .  .  Power  spectral  density, 
w 


<*> . Radial  frequency. 

. Fundamental  frequency. 


o  „  .  .  .  .  Cut-off  frequency. 


Discrete  frequency. 


CHAPTER  I 


.  INTRODUCTION 


1.1.  Objectives. 

^  -rl 

This- report  presents  the  results  of  a  mathematical  analysis' of 

V 

pressure,  flow  rate,  and  acoustic" power  levels  at  various  points  in  an 
idealized  piping  system  which  is  subjected  to  statistical  inputs.  Suc¬ 
cessful  adaptation  and  use  of  the  results  will  aid  in  the  prevention  of 
excessive  liquidborne  noise  in  shipboard  piping  systems. 


1.2.  Method  of  Attack. 

After  careful  review  of  existing  literature  and  technical  reports, 
an  analysis  was  begun  to  determine  the  feasibility  of  predicting  the 
response  v.f  piping  systems  to  liquidborne  noise.  The  analysis  was 
analytical  in  nature  and  used  the  results  of  previous  scientific  inves¬ 
tigations.  At  various  points  in  the  analysis  certain  assumptions  are 
made  on  the  basis  of  engineering  judgment.  These  points  are  clearly 
stated  and  discussed  in  this  report. 

It  is  shown  that  the  pressure,  flow  rate,  and  power  at  any  point 
in  a  complex  piping  system  can  be  found  in  terms  of  the  same  quantities 
at  some  more  convenient  location  in  the  piping  system.  A  digital  com¬ 
puter  is  needed  for  the  complex  systems.  The  results  are  expressible 
in  terms  of  continuous  spectrums  (i.e.  spectral  densities)  or  in  terms 
of  discrete  spectrums  and  bandwidths  as  are  obtained  from  wave  analyzers. 


1 


1.3.  Content  of  This  Report. 


This  report  contains  background  information  about  the  physical 
system  with  the  idealizing  assumptions  and  the  resulting  equations 
which  describe  the  response.  Further,  background  information  on 
mathematical  statistics  is  included  to  clarify  terminology.  The  mathe¬ 
matical  symbols  and  terminology  used  is  that  acceptable  as  "American 
Standards"^  when  applicable. 

The  analysis  follows  in  which  the  physical  system. is  analyzed 
*  ; 

using  the  mathematical  techniques*  The  results  will  be  discussed  else¬ 
where  in  this  report. 


*Numbers  in  parentheses  refer  to  references  in  the  Bibliography. 
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i  r°+T 

?T  /  X(t)  dt  , 

J  t~-T 


i  rto+T 

'  T  Jt»-T  * 


(t)  cos  nwjt  dt  , 


(2.5) 


,  .  2  f  0 


t«+T 


x(t)  sin  nojjt  dt 


*  .. 


By  further  manipulation  the  equations  (2.  3)  and  (2.  5)  can  be  placed 
in  "exponential  form"  as  follows: 


jnw.t 


V*  Jnwi 

X(t)  =  Z  Cne 


(2.6) 


n=-co 


where 


f  tn+ T 

1  .  i  I  0  -jnw.t 

'n  =  2  K  "  Jbn>  =  2T  /  *(*)  e  dt  *  (2. 

•'tft-T 


:n  “  5T  X  pT^l  - 

*-  s=jnu1 


(2.8) 


C0  =  a0  ’ 


(2.9) 


The  function  x^,(t)  is  defined  as 

xT(t)  -  x(t)  ,  t0-T<t<t0+T  ; 


Xfj«(t)  =  0  at  all  other  points. 


(2. 10) 


The  Laplace  operation  in  Equation  (2.8)  is  the  "one-sided" 
Laplace  operation. 


2.3.  Transfer  Functions. 


A  system  represented  by  a  linear  differential  equation 

L  [y(tj]  =  x(t)  '  (2.11) 

can  be  transformed  and  written  as 

Y(s)  3  H(s)  X(s)  .  (2.12) 

- 

The  solution  of  (2.  12)  is  ' 


y(t) 


s  JL  [Y(S)]  =  I  h(T)  X(t  -  T)  dT  . 

J  -00 


(2.  13) 

m 

Thus,  if  the' input  x(t)  is  known,  the  output  y(t)  is  found  by  an  inte¬ 
gral  operation  on  the  input  using  the  proper  weighting  system  function 
h(r)  . 

The  functipn  H(s)  is  the  transfer  function  and  in  some  cases  may 
be  an  impedence. 


2.4.  Output  Responses. 

When  x(t)  is  a  prescribed  input,  the  various  output  responses 
are  obtained. 

(a)  Impulse  response;  x(t)  =  6(t) 

y(t)  ■  h(t) 

Y(s)  =  H(s) 

(b)  Step  response;  x(t)  3  u(t) 


y(t) 


f* 

=  |  h(r)  dr 

Jo 


Y(s)  3  j  H(s) 


(2.  14) 


(2.  15) 
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(c)  Frequency  response;  x(t)  »  sin  ut 


JQ 


h(T)  sin  w(t  -  t)  dr 


Y<s>  58  tt-*  H<s> 

s  +  w 

where  by  common  acceptance 

H(jcj)  -•  frequency  response. 


(2. 16) 


(2.  17) 


It  is  assumed  that  the  effect  on  y(t)  due  to  the  poles  of  H(s).  has 
decayed  as  a  transient  and  that  only  s  =  ±  ju  will  be  effective. 


2.5.  Statistical  Concepts  and  Correlation. 

The  correlation  between  two  continuous  processes  or  motions 
is(3,4,5,6) 


012(t)  =  Hrn^  yj,  j_T  xx(t)  x2(t  +  t)  dt  =  x^t)  x2(t  +  t)  ^  ^ 

in  which  the  process  is  stationary  and  ergodic.  The  above  operation 
shown  is  denoted  by  Xj(t)  x2(t  +  t)  and  is  called  the  "time"  average. 
The  time  average  is  the  scan  of  Xj(t)  by  another  function  x2(t)  by 
shifting  an  amount  t  and  averaging  over  a  sufficiently  large  period 
2T  .  012(t)  is  the  cross-correlation  between  x^t)  and  x2(t)  . 

When  Xj  as  x2  ,  Equation  (2.  18)  yields  the  autocorrelation 


function 


r  x,(t)xl(t 


0n(T)  =  j_T  XjfDXjCt 

The  mean  square  value  is  when  t  =  0  ,  or 


+  t)  dt  . 


(2.  19) 


_  *  r  T 

t>U<0>  *  xl«  V£^3T  J.t 


Lm 


(2.20) 
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The  autocorrelation  ie  an  even  function  in  r  . 

Correlation  functions  are  useful  in  any  general  approach  to 
analysis  of  aperiodic  processes  and  are  instrumental  in  the  frequency 
description  of  a  random  process. 

2.  6.  Spectra  and  Spectral  Densities. 

The  spectrum  of  a  function  is  a  plot  of  the  Fourier  coefficients 
cn  versus  the  frequency  .  For  periodic  functions  the  spectrum  is 
a  series  of  discrete  points  at  specific  frequencies,  and  for  appjiodic 
functions  the  spectrum  becomes  a  continuous  plot  as  the  period  be¬ 
comes  infinite^2* 4*  This  latter  statement  is  equivalent  to  saying 
that  the  increment  of  frequency  approaches  zero. 

The  spectrum  mostly  used  is  that  historically  associated  with 
power  calculations.  Power  calculation  results  are  expressed  as  the 
product  of  two  functions  as  in  Equations  (2. 18)  and  (2.  19)  in  which  the 
average  value  is  obtained.  The  spectrum  of  the  power  is  expressed  as 
a  function  of  the  square  of  the  Fourier  coefficients. 

A  random  vibration  can  be  considered  to  be  the  sum  of  a  large 
number  of  harmonic  vibrations  of  appropriate  amplitude  and  phase. 

The  total  power  is  the  sum  of  the  power  of  each  component  of  the 
vibration.  This  power  is  distributed  as  a  function  of  the  frequency; 
therefore,  the  amount  of  power  associated  with  each  frequency  is  of 
interest  and  is  called  the  power  spectral  density.  The  accepted 
mathematical  definition  of  power  spectral  density  for  an  ergodic  pro¬ 
cess  x(t)  ^  is 

_  r  T 

G(f)  =  lim  i  I  x(t)  est  dt  ,  (2.21) 

T-oo  1  J-  T 
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where 


s  =  2irjf,  w  >  0  in  (2.21),  and’  f  =  ^  . 

Applying  Equations  (2.  3)  through  (2.  10)  to  Equation  (2.21)  with 


w  =  nwj  ,  f  =  nfj  ,  and  ^ 


(2.22) 


there  is  obtained 


G(f)  =  lim 
T-n 


x(t)  cos  nw 


(t)  sin  nwj  t  dt 


which  becomes 


G(fj  =  lim  i  |Tan  +  jTb  2  =  lim  I  c  f  .  (2.23) 
T-*oo  1  I  n  n  T— oo  *1  I  nl 

Thus,  as  a  function  becomes  random,  the  power  spectral  density  will 

be  given  by  Equation  (2.23). 

A  useful  relationship  for  statistical  signals  is  obtained  from  the 


equality^ 


CT-T 

J-T  X‘ 


(t)  x2(t  +  r)  dt  - 


/oo 

•oo 


xT1(t)  xT2(t  +  r)  dt  .  (2.24) 


After  substitution  of  the  two-sided  Laplace  integral  for  x^g (t  +  t) 
and  the  inverse  Laplace  integral  X^fs)  in  the  right  side,  dividing 
both  sides  by  2T  ,  taking  the  time  average  of  both  sides,  integrating 
the  left  side,  and  letting  T  go  to  infinity.  Equation  (2.  24)  becomes 

T  fal2^T^1  ~  $12^s^  *  "2T"  ^Tl^”8^  ^T2^  '  (2.25) 

L  J  T-*oo 


And  if  x  j(t)  =  x2(t) 


£  TpnWl  “  fuC)  »  lim  JT  XTi(-*)XT1(«)  •  (2-26) 

L  T  **oo 


The  spectral  density  $n(j")  is  found  by  replacing  s  by  jw  in 
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Equation  (2.  26). 

The  power  spectral  density  as  given  by  Equation  (2.  21)  is 
equivalent  to  • 

W)  c  lim  m  XT(s)XT(-s)  *  lim  ~|xT(s)|2  (2.27) 

T-oo  1  1  1  T-oo  1,1  1 

*  4 

Where  s  =  2irjf  .  Corftparison  of  (2.  26)  and72.  27)  shows  that  the 
power  spectral  density  is 

G(fj  «  2$u(f)  .  .  (2.28) 

Modern  control  systems  textbooks^’  ^  define  the  spectral 

density*  as  $j^(jw)  while  the  American  Standards ^  use  G(f)  .  No 
particular  difficulty  occurs  in  this  if  one  is  consistent,  because  most 
re*sults  will  be  normalized  in  the  same  manner.  The  factor  of  2  is 
accounted  for  in  that  $n(f)  allows  negative  f  while  G(f)  does  not. 
Thus,  both  quantities  when  integrated  over  the  frequency  spectrum 
must  give  the  same  average  power. 

Equation  (2.  23)  gives  the  relation  between  discrete  and  continu¬ 
ous  signals.  Using  Equation  (2.  28), 

$u(f)  »  Urn,  |  cn|  *  lim  5^- 1  cj  .  (2.29) 

Equations  (2.  23)  and  (2.  29)  are  more  symbolic  than  useful  and  show 
how  the  spectral  density  may  be  found  as  a  limiting  procedure  when 
the  frequency  fj  —  0  . 

2.7.'  Output  Response  for  Statistical  Inputs. 

A  system  with  multiple  random  inputs  is  shown  in  Fig.  1.  The 


*  Since  any  function  can  be  resolved  into  a  spectrum,  the  historic 
power  connotation  has  been  dropped  and  the  term  spectral  density  is 
used  for  any  function.  Thus,  force  may  have  a  spectral  density. 
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Fig.  1.  System  with  Multiple  Random  Inputs 


inputs  are  denoted  by  Xj(t)  and  the  respective  system  transfer  func¬ 
tions  by  Hj(s)  ,  i  =  1,  2,  3,  •••  ,  n  .  After  truncating  the  inputs  as 
in  the  previous  section,  the  following  equations  can  be  written: 

n 

XT0<S>  =  Z  Hi<s>XTi<s>  :  P-30) 

i=l 


and 


n 


XT0^"S^  =  2  XTk^"8)  • 


(2.31) 


k-1 


Multiplying  Equations  (2.  30)  and  (2.  31)  and  dividing  the  results 
by  2T  , 


n  n 


xto(s)  XT0^~S^  s  ~TF  Hi(s)  Hk(~s)  XTi(s)  XTk(-s) 

(2.  32) 


If  the  time  average  is  taken  over  the  above  equation  and  T  is  allowed 
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to  approach  infinity, 

n  n 

$0o(j“)  =  2  2  HkH")  $ki(j“)  (2.33) 

i=l  k=l 

where  s  has  been  replaced  by  ja>  . 

#  • 

Equation  (2.  33)  can  be  used  to  obtain  the"spectral  density  of  the 
output  in  terms  of  the  spectral  and  cross -spectral  densities  of  the  in¬ 
puts  and  the  system  functions.  As  an  illustrative  example,  a  single 

input  system  is  analyzed. 

» 

Single  input: 

1  1 

12  H^ju)  Hk(-ju>)  $ k.(j«) 
i=l  k=l 

1 

-  2  Hi0“>  Hi<-M  $UG“) 

■  i-  1 

>  Hj(j«)  Hj(-M  f  n(j«) 

-  |HjU«)|2  $  n(M  •  (2.34) 


2.8.  Spectral  Density  and  Amplitude- Frequency. 


An  amplitude -frequency  plot  of  a  process  as  obtained  from  a 


wave  analyzer  gives  the  absolute  value  of  the  Fourier  coefficients  cR 
for  a  bandwidth  Aw  with  center  frequency  u>n  . 

From  Equation  (2.  26),  for  r  =  0  and  s  =  jw  , 


(2.35) 


If  Xjft)  is  represented  by  the  complex  Fourier  series, 


x,(t)‘ 


lcnl 

— T 


(2.36) 


Then,  for  the  bandwidth  Aw  with  center  frequency  wr  ,  Equation 


(2.  35)  becomes 

V2 

T 


i /'■;>, .ov 

W_ - rr- 


dw  . 


(2.37) 


n  ~T 

Carrying  out  the  integration, 

1 2 

IT  IC  i 

,  =  I  n  I 
'll^V 


•  •  • 


•/* .  j '»/ 


t  ie  r 

*i,<K>  =  -7T 


(2.  38) 


Thus,  Equation  (2.  34)  of  the  preceding  section  can  be  written 

No  ■  lHi<M  Ni  <2-39> 

where  |cn|j  is  the  amplitude  of  the  input  Xj(t)  for  some  bandwidth 
Aw  with  center  frequency  wr  ;  |  Hj(jwn)|  is  the  absolute  value  of  the 
transfer  function  evaluated  at  w  **  wR  ;  and  |cn|0  is  the  amplitude  of 
the  output  xQ(t)  for  the  bandwidth  Aw  with  center  frequency  wr  . 


The  Physical  System 


2.9.  Introduction. 

The  following  material  contains  the  derivation  of  the  differential 
equations  describing  compressible,  turbulent  flow  in  a  non-rigid  cylin¬ 
der.  The  derivation  is  based  on  certain  idealizing  conditions  that  are 
pointed  out  and  later  discussed.  Using  Laplace  transform  methods, 
the  frequency  response  solutions  of  the  describing  equations  are  ob¬ 
tained.  The  results  of  this  investigation  are  applied  in  the  next  chap¬ 
ter  to  determine  the  response  of  a  complex  piping  system,  such  as 
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shown  in  Fig,  2, to  random  inputs. 


2,  10,  Derivation  of  Describing  Equations. 

Fig.  3  depicts  a  liquid  filled  cylinder  with  non-rigid  walls  con¬ 
taining  compressible,  turbulent  flow.  The  necessary  idealizing  con¬ 
ditions  follow  and  hereafter  shall  be  referred  to  by  number: 

1.  The  flow  is  one -dimensional. 

2.  The  frictional  resistance  is  a  function  of  some  nonlinear 

m 

operation  on  the  fluid  velocity  and  system  parameteNS'  It 

>  *» 

can  be  expressed  as  the  product  of  a  constant  and  the  instan¬ 
taneous  flow  rate  variation. 

3.  The  energy  dissipation  is  due  to  heat  transfer  and  is 
accounted  for  by  Condition  2. 

4.  The  rate  of  change  in  fluid  velocity  with  respect  to  length  is 
negligible  as  compared  to  its  rate  of  change  with  respect  to 
time. 

5.  The  instantaneous  pressure  is  constant  over  the  cross  section. 

6.  The  rate  of  change  in  mass  density  with  respect  to  length  is 
negligible  as  compared  to  its  rate  of  change  with  respect  to 
time. 

7.  The  velocity  of  wave  propagation  in  the  liquid  cylinder  is 
constant. 

Using  Conditions  1,  2,  3,  4,  and  6,  the  fluid  dynamical  relation¬ 
ships  completely  describing  the  system  are^*^: 

Equation  of  Motion 


&ut  ,  9p. 
"ST  '  p  Tic  + 


(2.  39) 
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terminating 

end 


..  ...  > 


Fig.  3.  The  Physical  System 


Equation  of  Continuity 


1  +  =  o  • 

r  QV  U  » 


T  at 

Equation  of  State 

*t  =  K'f 


(2.40) 


(2.41) 


where  K'  is  the  bulk  modulus  of  the  system,  i.  e. ,  the  fluid  and  cylin¬ 
der  walls  combined. 

Equations  (2.  40)  and  (2.41)  are  combined  to  give 


.  i 

‘  a*  T  at 


i=o. 


(2.42) 


Integrating  Equations  (2.  39)  and  (2.  42)  over  a  control  volume  A6x  , 


®Pt  P  8qt  n 

“"ST  +  &”8T  +  Bqt  s  0  * 


(2.43) 


and  ' 


»q+  a  8pt 

-TET  +  r-T5T’  0  • 

When  the  volume  flow  rate  is  constant  at  a  cross  section. 


(2.44) 


q^  *  q  ,  Equation  (2.  43)  reduces  to  the  slope  of  the  pressure  grade 
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line,  or 


|  =  Bq n  = 

“  2DA* 


(2.45) 


Then, 


.  — 2-n 

b  =  iea-y- 

2DA^ 


Pf 


iqn 


(2.46) 


where  p^  is  the  pressure  needed  to  overcome  the  frictional  resistance 

of  the  pipe  as  determined  from  elementary  fluid  mechanics  and  q  is 

the  mean  flow.  The  value  of  n  is  estimated  from  experience  or 

* 

experiment  and  should  range  from  1.  65  to  2.  05. 

The  nonlinear  term  of  Equation  (2.  43)  is  rewritten  and  expanded 
as  follows: 

q?  *  (q  +  q)n 
•  qn(!  +  f-)n 

»  q"  Q  +  "(4-)  +  n*n2;,  *>  (A)2  +  •••]  .  (2.47) 

The  ratio  -S—  is  small,  and  the  series  in  Equation  (2. 47)  is  suffi¬ 
ciently  approximated  by  its  first  two  terms.  Then,  Equation  (2.  43) 
becomes 


9pt 

"5F  T  r"5T 


pf 


npf 

_L  +  - -  q  =  0  . 

iq 


(2.48) 


Since  qt  =  q  +  q  and  pt  =  p  +  p  ,  Equations  (2.  44)  and  (2.  48) 
can  be  reduced  to  the  describing  differential  equations  for  the  varia- 
tions  in  flow  rate  and  pressure.  They  are: 


.|H+P^+^q  =  0; 


ax 


X  at 


*q 


(2.49) 


and 


- 1  1 "  ■  1 « ■.RTiir-wrc1*1 
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+ 


A  9p 

RTf  ' 


0  . 


(2.50) 


It  is  noted  that  the  above  equations  can  be  combined  to  form  the  classi¬ 
cal  acoustical  wave  equation  including  a  dissipation  term. 


2.11.  Idealizing  Conditions. 

i  ’ 

The  idealizing  conditions  listed  in  the  previous  section  are  dis¬ 
cussed  and  validated  to  varying  degrees  as  follows: 

1.  In  any  real  piping  system,  a  propagated  wave  is  never  one- 

dimensional*  because  the  elasticity  of  the  pipe  walls  allows  radial 

% 

motion.  But  over  a  certain  frequency  range  from  0  to  w  radians 
per  second  there  exists  only  one  propagated  wave,  or  one  mode  of 
propagation,  which  is  the  axial,  or  plane  wave,  mode.  This  mode  is 
denoted  the  (0,0)  mode.  The  cut-off  frequency  u  at  which  the 
first  radial  mode  of  propagation  occurs  £the  (0,  1)  modej  for  non- 
dissipative  conditions  can  be  determined  by  the  following  equations 
given  by  Jacobi^1 


JQ(Z>  _  2p2b  2Eb 

zJl<z>  “  ^  '  PlDc2z2  1 


(2.51; 


and 


-  2c 

wc  “  UZ01 


(2.52) 


where  z  is  some  characteristic  value,  Jq  and  J,  are  Bessel 
functions  of  the  first  kind,  and  Pj  and  pg  are  the  mass  densities  of 
the  fluid  and  pipe  material,  respectively.  The  first  real  characteristic 
value  Zqi  that  satisfies  Equation  (2.  51)  is  substituted  into  Equation 
(2.  52)  to  obtain  w  .  For  8"  standard  steel  pipe,  f  a  6900  cycles 
.  per  second.  For  smaller  sizes  of  pipe  the  value  of  w  increases  sub- 
stantially.  To  the  best  knowledge  of  the  author,  the  subject  of  cut-off 


*.ip»  . 
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frequencies  for  dissipative  conditions  is  nonexistent  in  the  literature. 

The  cut-off  frequency  should  be  somewhat  higher  for  dissipative  condi¬ 
tions  than  that  given  by  Equations  (2.51)  and  (2.  52).  It  is  concluded 
that  for  a  relatively  rigid  pipe,  the  one-dimensional  condition  is  valid 

for  frequencies  up  to  the  cut-off  frequency. 

.  *+ 

2.  The  perturbation  process  used  in  the  derivation  of  the  des¬ 
cribing  equation  to  linearize  the  frictional  dissipation  term  was  proposed 
and  verified  experimentally  for  low  frequencies  by  Waller^1^.  There 
seems  to  b&no  good  reason  for  discounting  its  validity  in  the  frequency 
range  of  one -dimensional  propagation. 

3.  The  effect  of  heat  transfer  on  a  propagated  wave  could  be 
considered  resulting  in  a  third  describing  differential  equation  as  done 
by  Browi/13),  but  it  is  believed  that  the  rather  empirical  fricitional 
dissipation  term  used  in  the  describing  equations  will  account  for  this 
effect. 

4.  The  condition  that  the  rate  of  change  in  fluid  velocity  with 
respect  to  length  is  negligible  has  been  used  for  many  years  in  des¬ 
cribing  weve  propagation  in  liquid  cylinders.  Phillips in  showing 
the  acoustic  wave  equation  adequate  for  describing  turbulent  flow  in 
liquid-filled  steel  pipes  has  shown  this  term  to  be  negligible. 

5.  If  the  one -dimensional  wave  condition  is  valid,  it  follows  that 
the  instantaneous  pressure  is  essentially  uniform  over  the  cross  sec¬ 
tion. 

6.  As  in  Condition  4,  the  rate  of  change  in  mass  density  with 
respect  to  length  has  classically  been  considered  negligible. 

7.  Strictly  speaking,  the  velocity  of  wave  propagation  in  a  liquid 
cylinder  with  non-rigid  walls  is  not  constant  with  frequency.  With 
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increasing  frequency,  the  flexural  vibrations  of  the  .vails  tend  to  im¬ 
pede  the  propagation  causing  a  decrease  in  the  wave  velocity.  Mathe¬ 
matical  expressions  to  determine  the  velocity  of  wave  propagation  as  a 
function  of  frequency  are  given  in  the  paper  by  Jacobi^ For  a 
relatively  rigid  pipe,  such  as  a  standard  steel  pipe,  the  velocity  of 
wave  propagation  in  the  pipe,  a  ,  is  essentially  constant,  and 


where 

% 

K.  -  KbE 

KD  +  bE  * 


(2.53) 

(2.54) 


Rather  than  calculating  the  velocity  of  propagation  by  Equations  (2.  53) 
and  (2.  54),  it  can  be  determined  for  different  pipe  sizes  and  fluids 
from  nomographs  available  in  many  papers  and  books^16\ 


Equations  (2.  55)  and  (2.  56)  can  be  solved  by  Laplace  transform 
methods  for  solutions  in  the  frequency  domain. 

Transforming  Equations  (2.  55)  and  (2.  56)  with  the  Laplace  inte- 


gral 

F(s)  =  1 

f  f(t)  e‘st  dt  ,  ~ 

J 

0 

ap(x,s) 

9x 

+  (R  +  sL)  Q(x,  s)  =  0  , 

(2.  60) 

and 

> 

l 

..  .:,****’* 

8Q(x,s) 

ax 

+  s  CP  (x,  s)  =  0 

(2.61) 

where 

p(x,  0)  and  q(x,  0)  are  zero.  Transforming  the  above  equa- 

tions  with  the  Laplace  integral 

- 

£<X>  -  1 

r  00 

1  f(x)  e  dx  , 

• 

J 

0 

-  X  f*(X,s) 

+  (R  +  sL)  §  (X,  s)  =  -  P(0.s)  , 

(2.  62) 

and 

-  X  Q(X,  s) 

+  sC  £  (X,  s)  -  -  Q(0,  s)  . 

(2.63) 

Solving  Equations  (2.  62)  and  (2.63)  simultaneously  for  £(X,s) 
and  Q(X,  s)  , 


f’(X.s)  =  Z  Q(r,s)(  ,  T.  )  +  P(r, s)(  ,  X  ,)  , 

x  -  Y  X  -  y 

and 

•  ft(X.s)  =  +  Q(r,  s)(  ,  X  >) 

c  X  -  y  X  -  y 


(2.64) 

(2.65) 


where 


y2  =  sC  (R  +  sL) 


I 


(2.  66) 
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„2  _  R  +  sL 
c  sC 


(2.  67) 


Evaluating  Equations  (2.  64)  and  (2.  65)  with  the  inversion  integral 

/e+joo 

e-joo 


A  Xt 

F(X)  eAt  dX  , 


-  P(x,s)  a  P(r,s)coshYX  +  Z  Q(r,  s)  sinh  Xx  ,  (2.68) 


and 

Q(x,s)  «  Q(r,  s)  cosh  yx  +  — -1  sinh  yx  . *-  ^ -^(2 ; *69) 

>  c 

% 

Tlie  frequency  response  of  the  system  is  obtained  by  replacing  s  by 
jw  in  Equations  (2.  66)  through  (2.  69). 

The  propagation  coefficient,  Y(jw)  ,  can  be  expressed  in  the 
form  a-  +  j/3  by  finding  the  roots  of  Equation  (2.  66).  The  attenuation 
constant,  a  ,  is  given  by 

1 

a  =  [^(a/h2  +  w2L2  -  wL)j  ^  .  (2.70) 

The  phase  constant,  j3  ,  is  given  by 

1 

j3  =  ( JFTjqj  +  UL)J  .  (2. 71) 

The  characteristic  impedance,  Z  (jw)  ,  can  also  be  expressed 
in  complex  form  by  finding  the  roots  of  Equation  (2.  67).  It  is  given  by 

zc  -  JJJ  (0  -  jo)  .  (2.  72) 

By  the  same  investigations  as  for  Y(jw)  and  Z  (jw)  ,  it  can  be 
shown  that 

Y(-jw)  =  a  -  j0  ,  (2.73) 

and 

Zc(-j“)  •  scte  +  jo)  .  (2.74) 
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CHAPTER  HI 


TRANSFER  FUNCTION  ANALYSIS 


3.  1.  Mathematical  Statement  of  the  Problem. 

The  system  which  is  symbolically  depicted  in  Fig.  2  encompasses 

most  any  situation  expected  in  shipboard  piping  systems.  The  noise 
*  l 

* 

source  at  n  is  a  pump  (or  valve)  which  discharges  into  a  complex 
system  with  typical  changes  in  size  as  indicated  by  A  ,  with  looped 
lines  (B  to  C),  with  stub  lines  (D),  and  a  termination  at  r  . 

The  fundamental  problem  is  to  specify  the  allowable  noise  levels 
of  pressure,  flow,  and  power  at  n  or  at  some  intermediate  point 
which  will  not  exceed  allowable  levels  at  the  terminus  r  for  a  multi¬ 
tude  of  piping  configurations  n  to  r  .  It  is  expected  that  in  the  over¬ 
all  analysis  it  will  be  advantageous  to  express  the  noise  levels  in  terms 
of  known  levels  of  pressure  and  flow  at  either  end  n  or  r  .  This 
requires  determination  of  transfer  relationships  for  the  various  ca3es. 

The  transformed  pressure,  flow,  and  power  at  a  point  x^  in  the 
pipe  i  are  defined  as  follows: 

p(Xi.s)  «  Epk(Xj,  s)  Pk(s)  ,  (3.1) 

and 

P(Xj,s)  -  Zqk(Xj.s)Qk(a)  ;  (3.2) 

QtXj,  B)  "  Eqk(Xi>  B)  Qk(B)  .  (3.  3) 

and 


Q(x.,s)  «  Ypk(x.,s)Pk(s) 


(3.4) 
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W^s)  =  Fpk(x..s)  |Pk(s) 


(3.5) 


W(xi,s)  “  ^(Xj.s)  |Qk(s)|  ; 


(3.6) 


where  k  can  be  either  the  source  n  or  the  terminus  r  . 


From  the  discussion  preceding  Equation  (2.  38)  it  is  evident  that 

the  transfer  functions  E  ,  Y  ,  Z  ,  and  F  are  needed  to  relate  the 

various  spectra.  Thus,  the  primary  objective  of  this  chapter  is  tp 

determine  thpse  transfer  functions  and  to  place  them  in  the  best  form 

* 

for  analysis. 

Before  proceeding  to  a  discussion  of  transfer  functions,  the 
method  of  determination  of  power  will  be  indicated. 

The  sound  power  at  any  point  x  is  found^  as 


*WT)  “  J?0 

kh  T-*< 


Z.TT  J 


p(x,  t)  q(x,  t+r)  dt 


(3.7) 


where  the  average  power 


W  =  «5  (0) 

ave  vpq'  ' 


“  ^  L 


$pq(j«)  d(ju>)  . 


(3.8) 


The  power  associated  with  any  given  frequency  is  found  by  obser¬ 
vation  of  the  spectrum  of  the  power  either  in  the  form  of  a  continuous 


plot. 


$w(j<*>)  “  2T  Q(x»  jw)  P(x,-j«)  , 

T-eo 

or  in  terms  of  the  discrete  Fourier  coefficients, 

».<*„>  -  [-  'ilr-H  ^  •  • 


(3.9) 


(3.  10) 
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From  the  preceding  discussion,  it  is  evident  that  once  P(x, s) 

and  Q(x,  s)  are  known  the  power  is  obtained  directly. 

« 

3.2.  Transfer  Function  Determination. 

To  determine  the  overall  transfer  functions  it  is  first  necessary 
to  work  with  a  specific  pipe  in  which  the  fundamental  parameters  E  , 
L  ,  and  C  are  constant  with  respect  to  that  pipe.  The  relationships 
between  the  boundary  conditions  and  pressure  and  flow  at  any  point  in 

m 

the  subsystem  i  ,  as  shown  in  Fig!  4,  are  given  by  Equation#^,!!) 
and  (3.  12)(U).  \ 

P(xitB)  =  P^s)  cosh  YjXj  +  ZciQ.(s)  sinh  ,  (3.11) 

and 

Pi(s) 

Q(x^,  s)  =  Q.(s)  cosh  YjXj  +  —  sinh  .  (3.12) 

ci 

4 


Fig.  4.  Notation  Sketch  for  the  Pipe  i 
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The  calculations  proceed  from  some  known  termination  by  suc¬ 
cessive  application  of  Equations  (3.  11)  and  (3.  12)  with  the  notation 


P^s)  =  P^b) 
Qi(0,s)  «  Q|{s) 


for  x.  =  0 


P^s)  =  Pi+1(0,s)  *  Pi+1(s) 

Qi(V8>  =  Qi+i(°»s)  a  pi+i(e) 


(3.13) 


for  Xj  =  i . . 


Further, 


» *-> 


P(XJ,S) 

wp T  "  z(xi's) 


(3.  14) 


from  which 


Pi(0,s) 

=  Zi(°'8>  =  Zi<8>  • 


(3.  15) 


Another  term  which  will  be  useful  in  transfer  function  determination 


is  the  reflection  coefficient 
l"*  Z,  -  Z__. 


i  *  zd  _  |  r  i  a 
r+  z ,  1 1  ii e  - 


(3.  16) 


A  further  consideration  is  that  the  frequency  response  is  all  that  is 
needed  for  spectral  analysis;  therefore,  in  all  that  is  to  follow 


s  =  jw  . 

For  the  section  i  ,  the  transfer  relations  are: 


(3.  17) 


P(*iJ“)  a  QjCjw)  [Zj  cosh  YjXj  +  ZcisinhTix.J  ,  (3.18) 

and 

Q(xi#jw)  =  Q^jw)  [cosh  +  ^-sinhTjxJ  .  (3.19) 

In  an  alternate  form^1®*  20^ , 
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Pfyjw)  *  Q,:(jw)  Zci  e 


■T. 


(3.  20) 


Qfy.  jw)  «  Q.(jw)  e 


ix:  1  “  n 


-2 y *** 


-n 


(3.21) 


Proper  application  of  Equations  (3.  15),  (3.  18),  and  (3.  19)  yields 
one  method  of  calculating  transfer  functions.  This  will  be  termed-..* 
Method  A.  A  second  method  uses  Equations  (3.  15),  (3.  20),  and  (3.  21). 
This  is  Method  B.  Both  methods  have  merit,  and  either  would  be  suf¬ 
ficient  provided  adequate  computer  facilities  are  available. 


3.  3.  Method  A  -  Transfer  Function  Analysis^). 

This  method  uses  Equations  (3. 15),  (3.  18),  and  (3.  19)  succes¬ 
sively  in  a  complex  system.  The  system  is  shown  in  Fig.  5,  in  which 
a  typical  section  is  shown  in  Fig.  4. 


Fig.  5.  A  Series  Piping  System 
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The  following  notation  1b  employed: 
BjCju)  -  cosh  YjXj  ; 

BJ(ja>)  =  cosh  Yjij  ; 

TjCjw)  =  Zc.  sinh  y.x.  ; 

,  T^(jw)  =  Zc.  sinh  y.i.  ; 

MjO)  ■  -J—  sinh  7jX.  ; 

1  Zci  11 

MJ(ja>)  =  sinh  y.i. 

^ci  1  1 


(3.  22) 


(3.23) 


j(3-.24) 


where  i  refers  to  the  pipe  number,  i  «  1,  2,  3,  ••  •  ,  m  .  Then, 
from  Equations  (3.  18)  and  (3.  19)  , 

PfrjJw)  -  B^jwJP^jw)  +  T.(ju>)  Q.(jw)  ,  (3.25) 

and 

-  BjO)  Q.(jW)  +  MjOwJPjOw)  .  (3.26) 

Since 

pk0“>  *  Bk-l0“)Pk-l<j“)  +.Tk-l<J“>V10“)  (3-27) 

‘  )  . 

and 

Qk(J“)  =  Bk-!(^)  Qfc-iO)  +  M^.1(ju,)Pk.1(jo,)  ,  (3.28) 

successive  substitution  for  Pk(j«)  and  Qk(j«)  ,  k  »  i,  i- 1,  •  •  • ,  2  , 
in  Equations  (3.  25)  and  (3.  26)  will  yield  P(x^,  jw)  and  Q(x^,  jw)  in  the 
form 

P(Xj,  j«)  =  G.(jo>)  Qr(j«)  , 

Gj(j«) 

P(*i.J«)  -  — z —  Pr(J«)  . 

r 
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Q(*i*  j«)  s  H.(jw)  Qr(jw)  , 

HjCJa) 

Q(x.,jw)  =  — Pr(jo>)  . 
r 


(3.  30) 


The  transfer  functions  G^(jw)  and  H<(jw)  ,  i  =  2,  3,  •  •  • ,  m  , 
are  given  by 


where 


°i  ■  BiGi-l  +  Ti  !  • 


Hi  ■  Bi  Hi- 1  +  Mi  Gi-1 


Gi-1  *  Bi-lGi-2  +  Tf-lHi-2  • 


(3.31) 

0 


(3.32) 


IJi-l  =  Bi-!  «i-2  +  Mi-i  J 

etc.  The  transfer  functions  for  i  -  1  at  some  point  Xj  are 
Gi  ■  ^  +  Ti  •  •  1 


H1  ■  Bi  +  TT  ■ 

r 


(3.33) 


For  Xj  =  , 


Gi  ■  7Tr  +  Ti  • 


Hi  "  Bi  +  r;  • 


(3.34) 


Then,  using  the  results  of  Equations  (3.  34)  in  Equations  (3.  31), 


the  transfer  functions  for  i  =  2  at  some  point  x2  are 
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This  theory  can  be  extended  to  any  number  of  series  components,  but 
the  transfer  function  expressions  become  increasingly  complicated  °s 
the  number  of  components  increases. 


The  analysis  of  a  complex  piping  system  using  Method  A  would 
proceed  back  from  some  termination  with  successive  calculation  of 
the  transfer  functions  at  all  points  and  frequencies  of  interest  includ¬ 
ing  the  junctions.  Thus,  Equations  (3.  31)  would  be  used  in  the  calcu¬ 
lation  of  the  transfer  functions  rather  than  expressions  such  as 
Equations  (3.  35).  Because  of  this,  Method  A  provides  an  easy  way  of 
computing  the  many  transfer  functions  needed  to  completely  analyze  a 
series  piping  system.  For  any  preliminary  analysis,  such  as  deter¬ 
mining  maximums  and  minimums,  or  for  an  insight  into  the  overall 
performance  of  the  system.  Method  A  gives  undesirable  results. 
These  things  are  more  easily  obtained  from  the  results  of  Method  B. 
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3.4.  Method  B  -  Transfer  Function  Analysis. 

This  method  makes  use  of  Equations  (3.  15),  (3.20),  and  (3.21) 
which  are  applied  to  one  pipe  to  obtain  the  results  in  a  form  where  the 
maxima  and  minima  are  more  easily  identified.  It  may  also  reduce 
the  computation  time  and  make  the  analysis  more  amenable. 

Taking  the  absolute  values  of  Equations  (3.  20)  and  (3. 21), 


|p(xiM  ■  K<M|Zcl|.  e 


OjXj  1  +  A*  +  2Ai  cos  ^ 

MIT  -CK 


(3.37) 


|Q(xi#  ja>)|  =  |  e  1  1  jj 


fl  +  -  2  Aj  cos  "1 2 

<  ■  ■  i  i  ■  i _ ■■  .  > 


il  -2lilcosei 


(3.38) 


where 


■  ini 


-Vi 


(3.39) 


6^  is  defined  in  Equation  (3. 16),  and 
tf'i  =  -  2^x.  +  ei  . 


(3.40) 


For  xfc  «  ,  Equations  (3.  37)  and  (3.  38)  become 

Ip  Mull  -  lo  (lull  Iz  I  eVk  f1  +  Ak  +  2Ak  C0S  K 

|Pk+lH  '  |QkM  |Zck|e  L1  +|Ik|  2  '  2|Ikl  cos 


-  2  Ik  008 


<y 


(3.41) 


Qk+lM  “  KM  e 


Vk 


+  A’ 2  -  2A'  cos  M  "I  7 

cr-ciw 


(3.  42) 


where  and  are  defined  by  Equations  (3.  39)  and  (3.  40)  with 


i  =  k  and  x^  = 
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Then, 


|P(Xi,  jw)| 

■  ^ 

|Qj(j")  |  , 

(3.43) 

|Q(X.,  jw)| 

|  QjO")  |  . 

(3.44) 

lpk+l<H 

■  Kpk !  • 

(3.45) 

and 

i  * 

K+i<M 

=  K* 
qk 

I  Qk(ju)l 

(3.46) 

where  Kp^  ,  ,  K^k  ,  and  K^k  are  defined  by  Equations  (3*:  37^, 

(3.38),  (3.  41},  and  (3.  42),  respectively.  Successive  substitution  of 

Equation  (3.46)  for  Qk(ja>)  ,  k  =  i,  i-1,  •  •  • ,  1  ,  in  Equations  (3.43) 

and  (3.  44)  yields 

i 

KPi  n  ^  KPi 

k=i- 1 

|Qr0«)| 

(3.47) 

and 

1 

Kqi  IT  Kik  ■  KH  Kqi 

k=i- 1 

|Q(x±,  jo»)|  = 

|Qr(M| 

(3.48) 

where  Kpi  and  R 

are  identical  to  |  G.(ju)|  and  J  H^(ju>)  J  , 

respectively. 

To  calculate  the  \arious  lk's  ,  the  point  impedances. 

zk(>)  - 

pk<J“> 

e*k 

(3.49) 

Qk(jw) 

e  * 

must  be  determined.  Equations  such  as  (3.41)  and  (3.  42)  give  the 
absolute  value  of  .  Equations  such  as  (3.  20)  and  (3.  21)  give  the 
phase  angle  as 
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«*k  “  ^o(k-l)  +  tan 


-1  Ak-1  8tn'(,k-l 


-  tan 


-i  -Ak-isin'l'k-i 


1  cos  K- 1 


(3.  50) 


where  is  the  phase  angle  of  .  Reducing  Equation 

(3.  50), 

.  _  .  •  ■  i  2Ak-isin'('k-i 


<*k  "  «*c(k- 1)  +  Un 


1  -  A1* 

1  Ak-1 


(3.51) 


Calculations  start  at  the  terminus  r  with  the  necessary  values 
known.  The  relationships  derived  in  this  section  are  used  to  compute 
the  K  functions  for  any  point  in  the  system.  With  these,  the  desired 


spectrums  can  be  determined. 


3.  5;  Piping  Systems  with  Parallel  Components. 


rod 


VCIBli 


In  Fig.  6  a  piping  system  is  shown  containing  pipes  in  parallel. 
The  parallel  components  are  not  necessarily  dimensionally  or  mater¬ 
ially  identical.  The  system  to  the  right  and  left  of  the  parallel  pipes 
consists  of  either  a  single  pipe  or  a  series  of  pipes. 


Fig.  6.  Piping  System  with  Parallel  Components 
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From  Section  3.  3, 

Pi(3«)  15  G!..1(jco)Qr(j «) 


Qi(j«)  11  Hl-i(3w)  Qr{.iw) 


(3.  52) 


(3.  53) 


To  analyze  the  complete  system,  Pi+1(jw)  and  'Q^+1(jw)  must  be 


found  in  terms  of  a  transfer  function  and  Qr(jw) 
The  following  notation  is  employed: 


-  cosh  ; 

% 

Tk(i“>  “  zck  sinh  Vk  ! 


'(fr.M) 

(3.55) 


MkCp)  R  —  sinh  Yk*k 

ck 


(3.56) 


where  k  refers  to  the  parallel  pipe  number,  k  =  1,  2,  •  •  • ,  m  .  The 


relationships 


Pi(3<*>)  =  P^P)  =  P2(j«)  "  a  pm(3«)  . 

Qi(jw)  =  Q2(j«) •+  Q2(j«)  +  *•'  +  Qm(j«)  . 
Pi+1(3«)  =  PnO»)  c  P2i(3«)  =  •••  =  Pm,(j“>) 


(3.57) 

(3.58) 

(3.59) 


Qi+l(3«)  =  QVW  +  Q2.(jw)  +  .*••  +  Qmi(3w) 

are  valid  where  the  subscripts  1,  2,  •  •  • ,  m  and  1',  2', 


(3.60) 


refer  to  the  points  shown  in  Fig.  6.  Using  Equations  (3.  57)  and  (3.  59), 
■  Pi+10»)  *  B'(ju)  P.(jo)  +  T^(joj)  Qjfju)  ;  1 


Pi+lCp)  “  B^(j«)  P  M  +  T£(jw)  Q2(jw)  ; 


-  B^0»)  P.(jW)  +  T^(j»)  Qm(j«) 


(3.61) 


Also, 
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Q !•(]“)  15  BJCjtoJQjOa)  +  M^jwJPjOw)  ; 

Q2»Cj“)  =  ■*'  M2^w)  Pi(jw)  ; 

>  (3.62) 


Qmt(jw)  -  B^(J«)  QJj»)  +  M^(jo)  P?(ju)  .  J 

Equations  (3.58),  (3.  60),  (3.  61),  and  (3.  62)  can  be  manipulated  as 

follows  to  form  rn  +  1  independent  equations.  Substitute  Equations 

... 

(3.52)  and  (3r53)  for  all  P.(joc)'s  and  Q^(j«)'s  .  In  Equations  (3.61) 
subtract  each  equation  from  the  first  equation.  Add  Equations  (3.  62) 
and  substitute  Equation  (3.  60).  The  results  are: 

T'jQjM  *  TJ,Q2(j^)  -  GJ.jtBJ  -B')Qr(jo)  j 

TiQjO^)  *  T3Q3(»  k  GJ.jCBJ  -  B' )  Qr(ja,)  ; 


TJQjC>)  -  T^CJa)  «  -  B'1)Qr(j«)  ; 


TiQjCJwi  -  T*nQjn(i»)  -  -  BpQr(ju)  W3. 63) 

Q.+10»)  ~  BJQjCw)  -  BJQ2U«)  -  BJQ3l»  -  ••• 

"  W-  .  "  “  Gi-l<Mi 


+  MJ  +  3V1.\  +  •••  -5*  + 


+  M^)Qr(j->); 


Q1(j,*»)  +  Q2(jU))  +  QS(J»>  +  •••  +  Qk(j«)  +  ••• 
+  Qm(iw)  j  H-.1Qra<*>) . 
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G|  ,(B'  -  B‘ )  -  T'H*  .  1 

G.'(ju)  =  B'  G.!  .  +  T'  — ^ .  (3.  68) 

i  ii-ii  Tj  +  Tjj 

If  the  parallel  components  are  dimensionally  and  materially  identical, 


Gj(.«)  =  B'j  c;_1  +  -i-TlHJ..!  % 


% 

*°0, 


■  BJHj.j  +  BM|Gj.j 


(3.69) 


(3. 70) 


for  any  number  of  pipes.  * 

The  transfer  functions  for  parallel  pipes  will  have  to  be  deter¬ 
mined  as  just  shown  for  both  Methods  A  and  B. 


3.  6.  Summary  of  Methods. 


The  results  of  the  three  previous  sections  provide  adequate 
means  by  which  the  transfer  functions  of  Equations  (3.1)  through 
(3.  6)  can  be  determined.  For  a  series  piping  system,  the  results  of 
Methods  A  and  B  are  related  to  these  transfer  functions  as  follows: 


For  pre.ssure, 


Method  Method 

A  B 


P(x,,ju>)  G'^jw) 
Epr<xi-'ui|  T  (>r  "  Z  (.M 


,  P(*,.ju>  GjO^)  K. 

Vxi'H  =  *  jf- 

pn 

.  P(*,.j“>j  .  _ 

V!xi'j“''!  “  i  “  lGi(jM)l  -  Kpi 


and  * 


Vxi'3“>| 


P(xj.  W  _  oi(i“)|  _  ^ 


(3.71) 


(3.72) 


(3.73) 


(3.74) 
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For  flow. 


|VvH  ■ 


Q(V  j«) 

FJFT 


PMj«] 


WM 


.  .  |Ht(j<-)  |  = 


|p^r 


Eqr(xi' » 


j  Q(x±,  ju') 
!’4r(jw) 


lHi(M|=  Kqi 


Eqn<V » 


i  Kai 

=  tgsa  ’  * 

qn 


For  power. 


.(XjJ«)|  = 


|W*“>I  ■ 


Fqr(Xi-j“>l 


Epr<xi-»l  I  V(Xi-j">l 

Epn<xi- H  lYpn<xi'M 

Zqr(*i.j“>llEqr(lci'j“'l 


Pqn<xi>H  ■  I  V^-^II  VXi’H 


(3. 75) 


(3.  76) 


(3.  77) 


(3.  78) 


(3.79) 


(3.80) 


(3.81) 


(3.82) 


Iffalnirfl  MS 


WiMH 


'^■»iVr  r  Atmiiiht'*  •  i  innift/nma* 


CHAPTER  IV  . 


H()f 

RESUME  *£/> 


*00 


% 


4.  1.  Background  Material. 


»ie 


The  background  material  contained  in  Chapter  II  .comprfsggs.^ 
much  of  the  literature  survey  undertaken  before  the  analysis  began. 
Continuous  evaluation  of  existing  literature  as  concerned  with  the 
analysis  i3  in  progress,  tut  the  results  of  Chapter  II  are  considered 
reasonably  valid  and  sufficient  to  begin  the  analysis. 


4.2.  Analysis. 

The  result?;  of  Chapter  II  indicate  that  the  analysis  of  a  piping 
system  with  statistical  inputs  will  best  proceed  from  a  "transfer 
function"  type  of  investigation.  In  Chapter  III,  expressions  for  deter¬ 
mining  the  transfer  functions*  of  series  and  parallel  systems  are  ob¬ 
tained  by  manipulation  of  the  system  describing  equations.  The  trans¬ 
fer  functions  given  by  Equations  (3.  71)  through  (3.  82)  are  used  in  a 
relationship  such  as  Equation  (2.  39)  to  relate  pressure,  flow  rate, 
or  power- at  :;ome  point  in  a  series  system  to  pressure,  flow  rate, 
or  power  at  another  point.  The  method  of  analysis  for  parallel  com¬ 
ponents  gives  transfer  functions  which  are  easily  employed  in  the 
overall  analysis. 

The  analysis  of  a  piping  system  will  require  electronic  compu¬ 
tation.  For  a  complete  investigation,  these  computations  will  be 
extensive,  and  Method  A  will  probably  be  more  appropriate  due  to  its 
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relative  simplicity  and  possible  shorter  computing  time.  Method  B 
provides  relationships  that  will  allow  the  maximum  or  minimum 
response  of  a  system  to  be  more  clearly  defined  before  numerical 
investigation.  It  also  puts  the  transfer  functions  in  a  form  that 
separates  the  effect  on  the  response  of  each  component. 

4.3.  Future  Effort.  REPR0[*UClBLE 

The  analysis  of  a  piping  system  by  Method  A  using  a  digital 

computer  is  now  in  its  latter  stages.  The  same  system  will- aWo  be 

*•'  i 

completely  analyzed  by  Method  B.  * 

Preliminary  investigations  into  the  possibility  of  predicting 
upper  and  lower  bounds  of  system  response  from  the  results  of 
Method  B  have  given  favorable  results.  It  is  foreseen  that  an  analysis 
such  as  this  may  make  a  complete  analysis  unnecessary.  The  results 
of  this  type  of  investigation  will  be  compared  with  a  complete  solution 


for  validation. 
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